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1.  Introduction. 

Mapping  has  always  been  a  favorite  tool  of  geologists  so  that,  with 
the  increased  accessibility  of  computers,  it  is  scarcely  surprising  that 
geologists  should  now  be  very  interested  in  studying  spatial  distributions 
mathematically.  Their  aim  is  to  use  some  mathematical  model  to  smooth  or 
contour  a  net;  of  data  points,  or  to  interpolate  between  them  or  to  calculate 
some  average  over  the  area.  (Two  recent  general  references  with  substantial 
bibliographies  are  [1],  [2].) 

There  are  a  number  of  possible  methods,  each  with  advocates,  and 
some  controversy  has  arisen.  In  the  U.S.,  the  most  popular  methods  have 

a 

assumed  that 

value  at  data  point  =  value  of  deterministic  function  +  random  error 

(1) 

and  consequent  application  of  the  least  squares  techniques.  In  South  Africa 
and  France  [ 3] ,  the  preference  seems  to  be  for 

value  at  data  point  =  value  at  a  point  on  a  random  function  (2) 

and  application  of  moving  averages.  In  actual  fact,  these  two  models  are  not 
really  distinct,  .t  seems  useful,  without  attempting  to  say  what  should 
be  used,  to  try  to  explain  the  theoretical  background  of  these  methods, 
because  this  is  not  available  in  an  elementary  exposition.  At  the  outset 
it  should  be  said  that  this  area  requires  considerable  statistical  research; 
it  has  so  far  been  mainly  the  preserve  of  very  theoretical  workers.  What 
little  there  is  in  the  literature  by  way  of  applications  has  been  written 
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by  geophysicist  so  perhaps  the  earth  sciences  as  a  whole  will  help  to  create 
a  new  class  of  statistical  methods  (see  scm0  of  the  references  in  [ <1,) 

2 .  Conceptual  Problems . 

Let  us  restrict  discussion  to  the  two  dimensional  problem.  Thus 
suppose  we  have  a  quantity  y  measured  at  points  (x^>  Xg)  in  the  plane, 
distributed  over  some  region  R.  Geometrically,  j  have  points  (with 
coordinates  in  three  dimensions)  (x^,  Xg,  y).  Through  this  cloud  of  points 
we  desire  to  put  a  smooth  surface.  It  is  hoped  that  this  surface  will  show 
how  y  would  vary  with  position  (x^,  Xg)  if  there  were  no  local  or  small  scale  ■ 
variation  the  effect  of  large  scale  disturbances  is  uspposea  to  be  incor¬ 
porated  into  the  smooth  surface,  which  may  be  called  to  trend  surface. 

Thus 

observable  surface  =  trend  surface  ->  variation  surface  (3) 

The  separation  into  two  parts  is  made  on  the  vague  basis  of  "large"  and 
"small"  scale  variation.  Our  interest  may  be  in  either  part.  The  trend 
surface  seems  appropriate  for  interpolation  and  averages,  the  second  for 
the  detection  of  anomalies,  e.g.,  discreteore  bodies.  It  is  clear  that, 
unless  we  can  introduce  more  prior  knowledge  of  the  particular  phenomenon  under 
study,  or  have  a  bright  idea  about  the  mechanism  after  a  first  look  et  the 
data,  any  analysis  must  be  quite  empirical.  It  follows  that  one  cannot  say 
any  method  is  "best",. 

In  section  3»  will  discuss  some  specific  mechanisms  which  might 
account  renl  iotically  for  some  geological  distr. ’'utio-ig .  Here  we  will 
continue  with  a  more  general  annlyois  of  (3).  If  (3)  is  compared  with  (l),  we 
no! Ice  that  the  two  pnrta  are  assumed  to  he  of  a  different  nature.  The  trend 
Is  assumed  to  he  a  deterministic  fnnci ion  nni  the  variation  to  be  random  or 
stochastic,  If  (3)  ts  compared  with  (."*),  there  in  no  dissection  but 
"random"  appears  again, 
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We  therefore  discuss  first  how  randomness  might  enter.  In  the 
first  place,  y  will  usually  be  measured  with  error,  i.e.,  if  'e  site 
is  revisited,  sampled  and  the  sample  measured,  a  different  value  would  be 
obtained.  This  is  well  understood  and  so  can  be  ignored  for  the  moment. 
Suppose  now  that  our  interest  is  only  in  the  region  R  that  was  sampled 
at  a  finite  number  of  points  -  we  will  later  consider  the  case  where  region 
R  itself  has  been  chosen  because  it  is  typical  of  many  similar  area,  'i.e., 
our  interest  is  in  the  population  from  which  R  was  somehow  chosen.  There 
are  infinitely  many  hypothetical  surfaces  passing  near  a  finite  number  of 
points.  Some  will  be  more  credible  geologically  than  others.  It  is  a 
standard  tactic  in  science  to  deal  statistically  with  an  infinite  hypothetical 
populations.  Explicitly,  we  must  define  the  class  of  admissable  surfaces 
and  suppose  our  particular  case  is  a  randomly  chosen  member  of  the  class. 

Apart  then  from  measurement  errors,  it  is  in  this  way  that  "randomness" 
enters  when  attention  is  confined  to  one  area. 

The  population  of  admisssable  surfaces  may  be  defined  in  several 
ways,  corresponding  to  (l)  to  (2).  We  will  return  to  this  in  a  moment. 

If  the  region  R  itself  has  been  chosen  as  typical  of  many 
other  regions,  even  complete  knowledge  of  the  surface  over  H  would  be 
insufficient  for  conclusions  about  the  other  legions.  It  is  necessary  to 
define  again  a  class  of  admissable  surfaces  and  to  rege^d  the  surface  over 
R  as  a  randomly  drawn  member  of  the  class  if  we  are  to  use  statistical 
arguments. 

To  define  a  class  of  admissable  surfaces,  we  are  guided  by  scientific 
plausibility  and  mathematical  and  computational  convenience.  The  older 
method  is  that  of  (l)  and  the  commonest,  choices  for  the  deterministic  function 
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are  (i)  an  algebraic  polynomial  in  and  (ii)  a  trigonometric 
polynomial,  i.e,,  the  first  few  terms  of  a  two  dimensional  Fourier  series. 

It  is  usually  implicit  in  papers  on  this  basis  that  the  choice  of  function, 
the  highest  degree  terms  used,  etc.,  should  make  the  "random  error"  have  x 
zero  mean,  be  uncorrelated  and  have  constant  variance.  For  example,  the  most 
thorough  study  in  the  algebraic  polynomial  case  -  Fraser  Grant  [4]  -  makes 
this  quite  explicit.  His  methods  may  be  repeated  using  trigonometric  poly¬ 
nomials,  or  any  other  set  of  functions.  In  detail,  they  depend  strongly 
on  having  observations  on  a  regular  grid  , 

However,  the  model  (1)  method  does  not  have  this 
requirement  and  so  has  appealed  to  geologists  -  but  the  "error"  should  be 
studied  if  the  model  is  to  be  checked  at  all  with  reality. 

The  "errors"  referred  to  in  the  last  paragraph  should  be  thought 
of  as  the  heights  of  a  random  surface  e(x1?  x^)  above  the  observing  points  - 
because  we  might  have  used  other  observing  points.  If  P'  =  (x|,  Xg) 

are  any  two  observing  points,  the  above  assumptions  read 


E(e(x|,  xp)  =  0  ,  E(e*(x|,  x£))  =  a2, 

E(e(xJ,  X,!, )  e(x£,  x£))  =  0  . 


(4) 


Clearly  the  last  assumption  is  unreasonable  P'  and  P"  are  close  together 
and  reasonable  if  they  are  far  apart  -  for  this  is  the  way  we  expect  spatial 
correlation  to  behave.  We  are  supposing  that  e(x^,  Xg)  describes  small  scale 
disturbances  only.  It  would  be  preferable  if  we  could  do  an  analysis  which 
replaced  the  last  two  assumptions  of  (4)  by 


E(e(P' )  e(P") )  -  c(P' ,  P") 


(5) 
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where  the  covariance  function  c(P’,  P")  decreased  as  P"  moved  away  fromP' . 
If  we  assume  that 


c(P’,  P")  =  c1(P'  -  P")  (6) 

i.e.,  that  c(P’,  P")  depends  only  on  the  vector  (x^  -  x|,  x"  -  xj)  the 
error  function  is  said  to  be  a  homogeneous  random  function.  If  we  assume 
more,  namely  that 


c(P',  P")  =  c2(r) 

where 

r  =  J(x^  -  x^+(x-x^ 


(7) 


the  error  function  is  an  isotropic  random  function,  because  the  covariance 
now  depends  on  the  length  but  not  the  direction  of  the  vector  separating 
P’  and  P".  In  neither  case  does  the  average  depend  on  the  absolute  position 
of  P‘  and  P".  These  are  the  generalizations  to  two  dimensions  of  the  notion 
of  stationarity,  now  conmon  in  time  series  analysis  -  which  is  the  cme  dimen¬ 
sional  equivalent  of  our  problem. 

The  spectral  methods  of  time  series  analysis  are  nc  familiar, 
particularly  to  geologist.'  who  have  followed  computer  applications  of  the 
Kansas  group.  The  weakness  of  the  model  (l)  work  has  come  from  ignoring  the 
two  dimensional  time  series  aspect  of  the  error  model  implied  in  model  (l). 

Having  come  this  far,  one  sees  that  it  may  be  unnatural  to 
separate  the  observable  surface  into  a  deterministic  and  a  random  function. 

For  the  small  and  large  sca.1*  disturbance  notions  find  a  natural  definition 
in  rpectrel  terms  -  high  and.  lov  frequencies.  And  t-hei-e  is  no  more  an  unnaturml 
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division  between  the  two.  It  may  be  better  to  define  the  class  of  admissable 
functions  in  terms  of  averages  over  the  class,  to  assume  homogeneity  or 
isotropy  and  to  base  the  analysis  on  this  assumption.  We  cannot,  in  this 
short  paper,  explain  how  the  analysis  would  then  go,  except  to  note  that 
moving  a. -erases  appear  naturally.  It  is  not  too  different  from  the 
trigonometric  equivalent  of  Fraser  Grant's  method.  Thus  it  is  very  awkward 
without  regularly  spaced  data.  Also  the  general  weak  assumptions  made  imply 
that  one  should  have  a  considerable  number  of  data  points.  This  is  one  form 
of  model  (2).  We  hope  to  deal  with  methods  based  on  this  model  in  detail 
at  a  later  date.  The  reader  will  find  a  list  of  relevent  references  at  the  end 
of  the  paper  [6]. 

Matheron  [5]  advocates  another  version  of  model  (2).  He  dees  not 
assume  that  the  averages  are  independent  of  the  absolute  position  but  that 
the  averages  of  changes  are.  Explicitly,  he  considers  random  surfaces 
f(xr  x2)  such  that,  for  all  x2),  (t^,  hg) 

|  2{[f(x1  +  hy  x2  +  h2)  -  f(x1(  x2»c)  =  \(h) 


where 


(9) 


-  K 


+  K 


This  is  n^t  the  same  as  saying  chat 


xj  -  f(x1(  X?)  -  f > o ,  0) 


(9) 


is  an  isotropic  random  function  because  (?)  does  not  allow  the  calculation 
of  c(f,  p")  v'-cn  ?*  and  F"  are  distinct.  The  extensive  literature  [  S] 
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generated  by  Matheron  and  his  co-workers,  particularly  Serra  has  only 
recently  come  to  our  attention  At  the  time  of  writing  it  is  not  clear 
how  assumption  (8),  even  with  a  particular  choice  for  the  function  •y(h), 
is  made  to  support  a  method  of  analysis.  Matheron  however  endorse  the  use 
of  moving  average  methods  which  he  attributes  to  Krige.  In  fact  he  uses 
"Kriging"  to  describe  the  act  of  taking  moving  averages,  Matheron  claims  his 
methods  do  not  require  regularly  spaced  data  points. 

To  conclude  this  section,  we  return  to  measurement  errors.  If  they 
are  present, one  must  adr1  to  the  random  surfaces  discussed  above  a  random 
surface  which  is  isotropic  and  in  which  the  spatial  correlation  falls  off 
very  rapidly,  i.e.,  the  equivalent  to  the  white-noise  of  time  series  analysis. 
3.  Some  mechanistic  models. 

There  are  two  relatively  simple  ways  of  building  probabilistic 
mechanisms  that  may  have  application  to  trend  surface  work,  i.e.,  generate 
classes  of  random  functions  of  geological  pertinence.  The  key  paper,  with 
a  large  bibliography,  is  Whittle  [6].  The  tract  by  Matern  [7]  is  also  a  basic 
aourse. 

An  example  of  the  first  method  would  be  the  following  construction, 
Suppose  ore  bodies  are  initially  distributed  randomly  in  the  plane  with 
different  sices  ar.d  that,  after  a  giver,  time,  their  "influence"  is  felt 
at  neighboring  points.  For  example,  their  contents  may  diffuse  outwards 
from  their  initial  positions  which  we  might  approximate  by  points.  If  the 
medium  is  usmogenous,  c(v,r)  could  stand  for  the  concentration  a*  a  point, 
a  vector  r  away  from  a  body  of  mass  w.  If  masses  fell  at  positions 

then  a  measurement  at  x,  y  (x),  free  of  measurement  error,  will  yield 

y(x)  ■=  r  (w^  x  -  Ri )  (10) 
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since  for  the  ith  mass,  r  =  x  -  R  .  The  sum  in  (10 )  is  over  all  the  masses. 
If  the  masses  are  distributed  over  the  whole  plane  by  a  Poisson  process  so 
that  the  mean  number  in  area  ■  A  is  \  A  and  if  the  probability  of  a  mass 
lying  in  (w,  w,  +  dw)  is  f(w)  dw  independently  of  its  position,  the 
stochastic  properties  of  (10)  may  be  worked  out  easily.  In  particular,  the 
process  is  homogeneous.  If  ?(w,  r)  depends  only  on  the  length  of  r, 
it  is  also  isotropic.  If  it  is  relevent,  then  interesting  problems  arise 
because  now  the  date  should  enable  us  to  estimate  \  and  f(w). 

This  model  could  also  be  used  for  linear  reefs  -  one  would  then 
distribute  lines  rather  than  points.  Another  generalization  would  arise  if 
the  Poisson  distribution  is  inadequate  and  has  to  be  replaced. 

The  statistical  treatment  of  data  on  these  models  is  in  its 
infancy.  Models  of  type  (10)  are  called  ir.o'd  n.’  pverage  representat ions  . 

Another  class  of  models  is  suggested  by  classical  mathematical 
physics  ir.  which  the  mathematical  description  usually  has  the  form 


y  (x) 


y  (x)  = 


e  (x } 


x  inside  R 


x  on 


s.  ^ 


(ID 


wr. ere  C  is  the  boundary  of  R  and  where  L  is  a  linear  operator.  Ft 

,  o  p  ?  o 

example,  if  y  =  temperature  at  equilibrium,  I.  *  :  5  /<>x^.  In 

usual  prn;-:  iral  oases,  the  temperature  distribution  over  *he  boundary 

determines  *  be  in*crioi  ‘  fn- r^i a! ■’ ir- :  ribut  ion .  Thus 


v  ( x  )  = 


b(x,  x 1 )  .  <*’ ) 


■■  V  • 


(12) 


v,hi'T’p  b  ( x ,  *  ’  )  In  reiatpil  In  (lif  Hiiu-ii’n  HuiHton  |.p  t.||p  problem  an1!  tn 
i'nl  ermlnetl  uniquely  by  |.  mu!  II.  Thun  If  »(*')  in  a  rutnlnm  funcU.n,  on  U, 
y(x)  Ui  n  rttri'inm  f'uni'l  Ion  In  (hr  interim,  |i  will  lie  naan  t.bat  (l.1)  Is  a^nln 
R  nii  tv  Ini'  itvnrnw  rrprnneiitJiHi.il  Mho  (in).  |i  in  merely  nhlftlnrri  in  a  rlirforsiit 
way,  Kluirl  flow  In  ilnrirr  ibeil  hIiuI  truly  no  |(  pppiiih  1 1 1 1  m  approach  aiioulrl  be 
r*t  1  b  .mil  for  nonir  proli I miin  , 

A  rm  if  1 1  more  illffiriill  nitwit, Inn  may  be  rnSnvnnt  t.o  many  unlmiirnl 
prnbl omn .  Above  wn  have  been  lalltlna  about,  bounilniy  value  probl emu  In  an 
inot.rnpH  inmiium  but  wH.b  random  boundary  vsiunii,  Tbr  cnnv»H'flB  a t t.ual  inn  -- 
a  niucUuiii  wboue  proper!  lan  van,  rain1  m!y  ami  alrplr  bntimlary  conditions  --  in 
muob  hardor  l.o  formulate  and  deni  with. 
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fir. me  Onimenta  rn 

"lea  variables  reglnnaliaces  et,  leur  pa tlmat.lon" 

Mftnnr.ri  et.  Tie,  l'/>0,  by  0,  M&theron. 

There  have  been  many  pub]  leations  by  Matheron  and  hia  colleague*. 

The  central  reference  in  the  above  hook  of  30')  pages  mathematic* ,  Thin 
book  in  not,  easy  reading  and  hAn  tio  point  of  contact  with  the  usual  literature  * 

in  fart,,  it  h&n  no  references  at  all,  except  to  text  books  on  stochastic 

processes.  It  seems  worthwhile  therefore  to  attempt  r  survey  of  ita  contents. 

Since  it  devciops  a  great.  deal  of  formalism,  we  try  only  to  describe  in  our 
own  words  arid  notation  the  essential  problem  Matheron  tackles  and  his  methods  - 
also  we  give  references  to  key  sources  in  the  Fnglish  statistical  literature. 
This  may  understate  the  range  and  depth  of  the  book.  If  so,  we  hope  it  provokei 
someone  elBe  to  translate  his  methods  at  greater  length. 

Let  f(x)  be  a  function  of  x  *  (x^>  x^,,,.^),  a  point  in 
n  dimensional  KiK'lldean  space.  We  would  like  to  estimate 

Q"Jf(x)dx,  (1) 

where  dx  =  dx^.-.dx^,  given  some  j  ^formation  on  f(x).  To  fix  ideas,  f(x) 
might  be  the  density  of  a  metal  at  a  point  x.  in  n  =  3  space.  Then  Q  is 
the  total  quantity  of  metal;  the  Integral  makes  sense  because  f(x)  is  usually 
zero  outside  some  region  R.  One  other  example  is  important.  Suppose  kD(x) 
is  z  o  for  x  outside  R  and  unity  for  x  in  R.  Then 

V  =  J  kR(x )  dx  (2) 

is  the  volume  of  Mir  Tl.  Af/Hu.  /ri  von  i  nt<  n'wn  t 1 . i  mo  values  of 
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k^x),  we  might  wiah  to  estimate*  V. 

To  evaluate  (l)  when  ('(*)  «  (>  outside  a  r«»# t > ’l.  H,  ah’1  wh*M  Hip 
filin’.  t.lnnal  form  of  f  ( n  )  is  given  is  a  matter  of  UU.*«rftt  1-  n  ■  nuinf  i  S  >•»  I 
Integration  If  the  rum  of  rU)  In  >•<  ni>-l  l<  «*  e.l .  inn  divides  It  lot"  M 
subregions  R  with  volumes  V’  evaluates  f(s)  at.  *i  t  vp  I  '*»  |  (mini  x  j  io 
sat'li  Rj  and  forms 

II  *  r  T(k  )  V  .  < 

1-1 


Fancier  method*  approximate  f  locally  hv  simpler  I'miicI  U'ns.  but  the  llnal 
formula  will  have  a  form  like  (i).  'Him  *rrm  made  t»  a  Imply  'l-ll*  «hd  hooks 
on  numerical  integration  suggest.  ways  of  estimating  this  usUic  one  a  ki  1 1  < 

of  the  formula  for  f(x).  Thla  la  all  classical  mat, hemal  lee .  If  nothing  In 
known  about.  f(x)  except.  Its  values  at,  Xji....x^.  nothin,,  can  hs  said  about, 
the  error  Q-!3. 

Two  examplen  show  an  interesting  result  (Altkrn  (19  19).  v*  •  pe  (I'tliM) 
Kendall  (19**2)).  In  one  dimension  t.he  use  of  poU.tc  oijcnlly  sicced 
h  units  apart  reveals  tp-  fMlowii.t,  facts: 


00 

T 


-CO 


_1_ 

72n 


-  -  (o+nh)4* 

e  h 
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very  well  even  for  h  as  large  as  twq  while 


n(l-(fl+nh)  ) 


estimates  j 


dx 


n(U(x-of ) 


(5) 


much  more  poorly  for  similar  values  of  h  .  As  Kendall  showed,  the  reason  for 
the  difference  in  t.he  behavior  is  t  hat  the  charset  erist  io  fvmc.icn  (Fourier 
transform),  e 
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m  J  f(x)  '  y  / 
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n  Min*  m  |n)  h  hh  iinlilnnnil  nnUiralos  <>f  ij.  Also 


vm  (n(t-) )  **  Mn'’("))  | Pi ;n ) I'* 
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in 

'Minn*  Hm  niim  In  ('))  In  from  -*  to  ••  exrludtnn  m  -  u,  and 

*(t)  -  J*  •  lKt  f  ( x )  .lx  (10) 

mm 

la  Hm  Fourier  trarmf.'im  of  f(x),  Tim  error  In  fl(n)  will  bo  amaller  on  the 
iviriin  the  iiirallar  (•>)  la.  F"r  fixed  h,  thin  meana  | ^. ( t ) (  tending  to  aero 
faat.ar  *•  j  t]  *<  »,  The  behavior  of  (U )  and  (S)  can  nev  be  predicted.  If 
o  la  not  taken  aa  a  random  variable,  Kendall  a hewed  h> v  to  ,{et  aimilar  reaulta 
by  noticing  that  fj(fl)  and  the  *rn  r  4-h(n)  are  periodic  funotlona  of  8, 
the  period  belt,*  h.  and  repi<  Rentable  ar  K.nu  ler  flerion 


I 


Finally,  the  variance  of  fl(h)  would  also  be  known  if  we  know 

at 

hJv)  -  J  f(x)  f(x>y)  dx  (11) 

-n* 

tv»i*  every  y  -  heraune  then  w«.  could  evaluate  (0).  Mat, heron  calls  (11) 

the  covnflnprnm_of  f(x).  He  has  given  the  n-dtmenn lonal  version  of  the 
above  refiult.ii  the  only  cimt, pen  required  are  nut.nl  tonal.  He  alno  defines  the 
drml  -  vnr  1 '  <i’j,ntn 


y(y ) 


r,  J  |  f  (x+v )  -  f(x)t  dx 


(1?) 


Thun 


y ( v )  t,  J  (f  (xqy)  -  2  r(x-fy)  f(x)  *  i’“(x))  dx, 

(13) 


'  »'(0)  -  g( y) 


However  there  are  fiiictl-ns  for  which  (11)  is  infinite  but  for  which  (12)  is 
nr  t  -  t.hic  in  t  he  reanor  why  Matheron  prefers  (12).  We  will  return  later  to 
them*  questions  and  to  the  key  quention  of  choosing  a  form  for  g(y)  or  'y(y) 

In  practice. 

Sampling  methodn  for  reiving  mathematical  problems  (we  have  Just 
reen  ar.  example)  are  now  called  Monte  Mario  methods  and  are  very  well  known  - 
aee,  e.g.,  the  bo'k  by  Hammers  ley  and  Har.decom.be  (lf^li).  They  all  depend  on 
rtatictical  theory.  The  ab-'  ve  method  Is  by  no  means  '  ne  only,  or  the  best, 
way  of  evaluating  integrals  by  sampling. 
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Certain  geometrical  problems  are  oloaely  related  to  the  above. 

For  example  Kendall  ( L9UO)  considered  estimating  the  area  of  a  region  R  by 
counting  the  number  of  lattice  point*  which  fall  in  it  when  the  lattice  la 
poll tinned  at  rrndom.  The  theory  Jn  an  above  using  the  function  k^(x)  defined 
for  equation  (2),  Buch  problems  are  of  greet  mathematical,  interest  and  can 
be  traced  through  hooka  and  papers  on  integral  Geometry  or  Geometrical 
Probability  -  see  e.g.,  Kendall  and  Moran  (19*3),  Moran  (19*9).  Mathcron  culls 

K(h)  -  J  k(<(x)  kR(xi-h)  dx  (1*0 

the  geometric  covariogiam. 

So  far  we  have  here  regarded  f(x)  aa  an  ordinary  function  -  the 
only  stochautic  element  has  been  the  sampling  process.  But  workers  in  many 
fields  including  geology  have  considered  the  case  where  f(x)  Is  a  random 
function  with  certain  properties  -  the  reaaon  why  and  the  basic  references 
have  been  given  in  the  body  of  the  paper. 

In  this  caae  even  a  complete  knowledge  of  f(x)  so  that  (l)  may 
be  evaluated  ia  not  enough.  The  value  of  Q  obtained  Is  .lust  one  possible 
value  out  of  an  eumamble  whose  average,  K(q),  we  seek.  But 

E(Q)  -  J  E(f(x))  dx, 

■  j  ffi(x)  dx,  (15) 

where 

E(f(x))  •  m(x)  .  (l6) 

Since  var(Q)  «=  E(Q2)  -  E(Q)2,  we  need 

E(Q2)  -  E  jj  f(x)  -(x* )  dx  dx’  . 
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But  setting  x'  ■  x  *■  y,  we  have 

F.(Q2)  -  E  J((  (f(x)  f  (x*y ))  dx)  dy,  (17) 

and  using  (11), 

E(Q2)  -  J  F  g(y)  dy  .  <lR) 

This  should  be  compared  with  (0), 

E(f/  (0))  -  hT  g(nh)  .  (1°) 

•h' 

The  relationship  is  closer  than  one  might  have  thought.  In  this  approach 
we  usually  ask  whether  our  under**'  and! ug  of  'lie  problem  eng;  eats  n  toim  to 
assume  f  >r,  not  E(p(y))but, 

cov(f(x),  f(x+y))«  h(i'(x)  f(x+y)),  (^0) 

called  the  covariance  function,  I'.juivalcntly  one  night  try  to  spcci»y  a 
new  versi  i,  of  (12) 

E (\ (y ) )  -  \  E  { f (x-»-y )  -  f(x)  f  (-1) 

unieh,  if  known,  implie#  the  value  of  (20)  when  that  exists. 

To  estimate  E(Q),  giver,  only  f(x)  at  x^,...,x^  we  might  consider 
the  estimator 

F  =  x1  f(xJL)  +  ...  +  XR  f(*r) 


(22) 
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where  the  Xj'n  lire  to  be  chosen  to  minimize 

l’(R  -  Q )'  -  mean  square  error  of  R  .  (.'i) 

nut 

!’’(R  -  Q)?  ••  ff((f(x)  f(x'))  dx  dx' 

-  2  ^  ).i  J*  K(f (xi)  f (x))  dx 

+>  T  Kj  E(f(x1)f(xj))  •  W 

Hence 

2 

..2|  t;(f(x .)  f(x))  ch: 
o  A  ^  i 

♦  2  f(x1)  f(x.) 

-  o  (in  1 , . . . , n )  (25) 

la  &  oet  of  "normal"  equations  for  the  .  They  may  be  solved  if  t 

i.  n 

covariance  function  (t?0)  is  known.  Mntheron  calls  thi6  method  Xrlging. 

There  la  aome  literature  rn  how  best  to  ohcoae  the  points  x,  in  v 

x 

where  the  ran ’^n  function  i*  <J  reive<l  .««  tliAt 


f 


1 

N 


N 


r  r(x.) 

i 


"beat"  approxirateo 

n  J  f(x)  dx  . 
V 


■  A  9  - 


Mu  |,-u  (1'JV/I  ijri;  ••  remit.  ,hic  1 1  Tying  equal  ■  pacing  when  V  ifl  an  interval. 

!  a  1  eti  tun,  lln.|H«  m,  1  'uhrzv'kl  (]r/C>)  have  studied  the  planar  cane , 

fl"  f'nr  wt*  hav  ■  cunli’i  <  ,1  only  the  est  iiration  rf  the  integral  (1)  - 
thin  ttoeinfi  t.n  hr  the  <  nly  problem  attacked  in  the  book  although,  the 
m  t  1 1  n  1 I  nil  i if 

“  J  p(x)  f(x)  .'x  ,  (c6) 


w)iere  p  ( x )  in  Known,  in  considered,  if  p(x)  >»  6(x-x  ),  the  Dirac  delta 

0 

fundi  <n,  P  ■  f(xo).  Tcx.r,  in  jnrticular  (26)  raises  the  problem  of 
estimatin'*  f(*)  at  unobserved  point?,  'ibis  problem  is  considered  in  our  first 
reference  to  [•'atueren ' ?  works.  His  approach  seems  to  be  the  following. 

Suppose  f (• )  is  observed  at  y  , ...,x^  and  required  at  Write 


f(xo)  «  f 


°n  f(xn} 


(27) 


ar.d  choose 


1 


t  hat 


*  O 

2(f (xc)  -  f(xn)V  u  min isur. 


(26) 


Ih.is  least  squares  problem  may  be  solved  if  one  kr.v the  covariance  function 
1.29),  The  coefficients  are  clearly  dependent  cn  x  *  so  that  (.??)  is  a  kind  of 
roving  average.  Mather^r.  also  calls  this  hriging. 

All  the  av-ove  theory  car.  be  rewritten  with,  f (• )  vector  valued 
ar.d  iAthercn  does  this. 


Thus  all  the  above  theory  is  quite  simple  ar.d  net  novel.  Its 
application  is  seer,  always  to  in.-.lv**  he  choice  of  a  function  -  the  variogra:., 


deri-vari -*grem  .  c'vi;  : oi  ee.  It 
recommend,  and  why!  fhr  ild  arc  ,:.e 


difficult  to  see  y.  st  what  Mat heron  world 
vi '  1  int;  >  ■  accept  he  challenge  to  put 
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Matheron's  work  into  the  English  literature,  this  is  the  point  that  needs 
most  attention.  It  is  clear  that  he  favors  functions  which  are  isotropic. 
Even  giver  this  restriction,  the  estimation  of  the  above  functions  from 
data  is  non-trivial  and  this  does  not  seem  to  be  attempted. 
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The  methods  for  analysing  data  in  one  dimension  (usually  time  )  are  highly 
developed.  However,  in  several  dimensions,  most  applied  workers  only  use  one  class 
of  methods  —  that  in  which  the  data  is  assumed  to  be  the  sum  of  a  deterministic- 
trend  and  an  uncorrelated  random  error.  The  more  general  model  with  a  spatially 
correlated  error  has  been  neglected  in  most  fields  --  oceanography  is  a  conspicuous 
exception.  This  neglect  is  partly  due  to  the  lack  of  expositions  for  the  applied 
workers.  However,  the  manner  in  which  much  geological  data  is  now  collected  does 
make  this  application  difficult. 

The  aim  of  this  paper  is  to  explain  the  relevance  of  the  possible  models 
and  methods  and  to  indicate  their  data  requirements.  > .n  append  La  bar  been  added 
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